Renormalization Group Flow of Hexatic Membranes 
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We investigate hexatic membranes embedded in Euclidean _D- dimensional space using a re- 
parametrization invariant formulation combined with exact renormalization group (RG) equations. 
An XY-model coupled to a fluid membrane, when integrated out, induces long-range interac- 
tions between curvatures described by a Polyakov term in the effective action. We evaluate the 
contributions of this term to the running surface tension, bending and Gaussian rigidities in the ap- 
proximation of vanishing disclination (vortex) fugacity. We find a non-Gaussian fixed-point where 
the membrane is crinkled and has a non-trivial fractal dimension. 



A. Introduction 

Hexatic membranes are a very interesting example 
of physical system which undergoes a non-trivial phase 
transition induced by a combination of geometrical 
and topological interactions |1,]. This system is even 
more interesting since it is characterized by local re- 
parametrization invariance, and is thus an example of 
a theory characterized by local symmetries with a renor- 
malization group (RG) flow exhibiting a non-Gaussian 
fixed point in two dimensions. 

In the continuum, Hexatic membranes are described as 
a fluid membrane coupled to an 0(2)-symmetric scalar 
field. The 0(2)-global symmetry represents in the con- 
tinuum the Zg-symmetry characteristic of the orienta- 
tional order of triangular or hexagonal lattices (from this 
last it takes the name) that melt [ij; the 0(2)-symmetric 
scalar field plays thus the role of the orientational order 
parameter. Hexatic membranes are effectively described 
by combining the fluid membrane and the 0(2)-model 
actions: 



5[N,r]= I d':cV9{t^ + ^K' + ^R} 



(1) 



In ([T]) Ka is the hexatic stiffness, fi is the surface tension, 
K is the bending rigidity and R is the Gaussian rigidity 
[l|. The action ^ is invariant under re-parametrization 
of the embedding field r : M^ — > R-^; 5^1, = 9^r • dj^r 
is the induced metric, K^ is the square of the trace of 
the extrinsic curvature, while R is the intrinsic Gaussian 
curvature. The scalar field, tangent to the surface, is 
constrained to have unit length N • N = 1. 

In this paper we study the RG flow of the couplings 
appearing in the action ([T]) by computing their beta func- 
tions for arbitrary embedding dimension D. Differently 
from previous approaches to hexatic membranes [2-4], 
we use a geometrical formulation of the model which is 
explicitly covariant and we apply to it the effective aver- 
age action formalism Q , which has been adapted to fluid 
membranes in [6] . The phase diagram of fluid membranes 
is modifled by the inclusion of the hexatic stiffness, which 
induces, in the RG flow, a non-Gaussian IR fixed-point 
at which the surface is crinkled, i.e. it has a non-trivial 
finite fractal dimension [2J. 



The 0(2)-action in (IT|) can also be seen as the general- 
ization on a membrane of the continuum action describing 
the AF-model. As it is well know, the AF-model is sub- 
ject to the Kosterlitz-Thouless topological phase transi- 
tion Q, mediated by the unbinding of vortices. In the 
context of hexatic membranes, disclinations [ll| are the 
relevant topological excitations; these interact by long- 
range Coulomb interactions, and in this way they affect 
the phase behavior of the model. Interestingly, also the 
Gaussian curvature plays a similar role, since it is a source 
of frustration for the field N, which, when parallel trans- 
ported around a closed loop, is shifted proportionally to 
the total Gaussian curvature enclosed. This is described 
by long-range Coulomb interactions between Gaussian 
curvatures in the form of a Polyakov term J y/^R-^R Q, 
resulting from the integration of the 0(2)-field. 

Even if, in the effective average action formalism, it is 
possible to account for the effect of topological defects 
without explicitly summing over them [9\, in this paper 
we will discard the contribution of disclinations and focus 
on the effect that the long-range Coulomb interactions 
between Gaussian curvatures have on the phase diagram. 
Our main interest is to understand the effect that non- 
local invariants have on the RG flow of the couplings of 
the local ones. As we will see, it is the hexatic stiffness 
Ka, seen as the coupling of the Polyakov term, that has 
the effect of changing the phase structure of fluid mem- 
branes by modifying the beta functions of the bending 
rigidity. 



Induced long-range curvature— curvature 
interactions 



At every point of the membrane we can introduce 
zwei-bein orthonormal vector fields ea — e^daV such 
that e"e^gai3 = Sab- Covariant derivatives of the zwei- 
beins define the spin connection uJaab through the relation 
VaBa — uiaab^b, which dcscribcs how the zwei-beins ro- 
tates under parallel transport. We can now write the vec- 
tor field N in the zwei-bein basis N — cos 6ei+ sin 6 62 , 
defining in this way the angle variable 9, in terms of which 



the action ([T]) becomes: 



+ ^Jd^X^{da9 + Wa){d^9 + Lj''). (2) 

The angle 9 parametrizes the one dimensional sphere 
0(2) = S^. Since the field 9 appears quadratically, we 
can complete the square to perform the Gaussian in- 
tegral obtaining in this way an effective action for the 
membrane alone. This amounts in shifting the field 
6 ^ 9 — J -g-VaO;" of the amount given by the solution of 
the tree-level equation of motion A9 + Vaw" = 0, where 
A = — 5"^Vq,V/3 is the covariant Laplacian. After using 
the relation between the Gaussian curvature and the spin 
connection, R — ^eabV^w'', one finds that the action ([2]) 



becomes: 



S[9,r] 



V9'- 



= I (fx 



^Jd^x^da9d^9 



Ka 

8 



(3) 



The step between ^ and (|31) is analogous to what 
happens in Liouville theory, where the Gaussian action 
/ Y^ (i(/)A0 + 0i?) gives rise to a Polyakov term with 
both classical and quantum contributions. In the Liou- 
ville case one makes the shift —>■(/) — J -^ i? obtained by 
solving the tree-level equation of motion A(j> + (j>R = 0. 
The underlining idea is that the tree-level equations of 
motion directly give the correct change of variables to 
complete the square in a Gaussian integral. 

The 0(2)-field has the effect of inducing long-range 
interactions between Gaussian curvature; these are rep- 
resented in ([3]) by a term proportional to the Polyakov 
effective action. At this point, as done in |3], one can 
integrate out the field 9 to obtain the hexatic membrane 
effective action: 

S,Mr]-Jd'x^{f^o + ^K'+~^R} 

+^jd^x^gR^R, (4) 

with renormalized jjlq, kq (the bending rigidity is not 
renormalized kq = '*oo) and finitely renormalized hexatic 
stiffness: 



Ka.o — Ka, 



1 
12^ 



(5) 



One can obtain ^ using the effective average action for- 
malism along the lines of [Ifl], where it was shown how 
the Polyakov action arises as field modes are integrated 
out, one by one, from the UV to IR, and resulting in ^. 



C. Flow equation and beta functions 

The effective average action is a functional that de- 
pends on an infrared scale k and that interpolates 




Figure 1: Beta function of tlie coupling at for (from bottom) 
Ka = 1,1,2. 



smoothly between the bare, or microscopic, action for 
fc — )• oo and the full effective action for fc — >■ 0; when local 
symmetries are present it is constructed using the back- 
ground field method. The adaptation of the formalism 
to fluid membranes has already been done Q to which 
we refer for more details. For a complementary study of 
polymerized membranes using the effective average action 
approach see [11]. One can construct a functional Ffc [r] 
such that it is a re-parametrization invariant functional 
of the embedding field r of the membrane. 

The main advantage of working with the effective av- 
erage action formalism is that we dispose of an ex- 
act RG flow equation that describes its scale derivative 
(i^logfc): 



dtVkir] 



1^^_ dtRk{A) 



Ti'\r]+R,{A) 



(6) 



In ^ Tf, [r] represents the Hessian obtained by expand- 
ing the effective average as r ^ r + i/*n* to second order 
in the fluctuation fields j^*. Here ly^, with i = D — 2, are 
the fluctuations in the direction of the n' normals vec- 
tors and the function i?^ (A) is the cutoff kernel. The 
derivation of the flow equation ^ is also given [6|] . 

To proceed we need to make a truncation ansatz for the 
effective average action in order to project the RG flow 
to a treatable finite dimensional subspace. We choose the 
scale dependent version of (U): 



TfcM 



d^'x 



v?{ 



Mfc 



HtK'^ + It 



KA,k 



] 

Jd^x^R^R, 



(7) 



where we introduced the running surface tension /i^, the 
running bending Kk and Gaussian k^ rigidities. 

The hexatic stiffness Kam does not renormalize, i.e. its 
beta function is zero: 



dtKA,k = . 



(8) 



One way to see this, following the arguments of [2!|, is 
to consider the Polyakov term in the effective action Q 
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Figure 2: The flow in the (a^, flk) plane for D = 3 and Ka ~ 
1. The Gaussian and non-Gaussian fixed points are marked 
dots. 



as the result of the integration of N free scalar fields; 



then N = Kao = Ka. 



12ir 



is preserved by fluctua- 



tions. Another way to reach the same conclusion is to 
consider the 9 scalar field as a non-linear sigma model 
on 0{2) — S^. The beta function of KA,k is then pro- 
portional to the Gaussian curvature of S^ |l2| . which is 
identically zero since the manifold is one-dimensional. A 
third argument is related to the flow of the Polyakov ac- 
tion. In [101 it is shown that the RG flow induced by a 
minimally coupled scalar, as is 0, induces a Polyakov term 
in the effective average action only in the limit fc — > 0; 
thus for non-zero k, the flow of KA,k, as extracted form 
the Polyakov term, is zero. In the following we will con- 
sider KA,k = Ka as a parameter. As in the standard 
Kosterlitz-Thouless topological phase transition, the flow 
of the hexatic stiffness is induced by topological defects 
(here the disclinations) [Ijl- We will discuss this in a 
future work |l4| . 

To obtain the projected flow for the couplings in ([7]) 
we proceed as in p| with the difference that we need to 
determine the contribution of the hexatic stiffness to the 
running couplings, i.e. we need to compute the terms 
proportional to / y/g, J ^/gK^ and / ^/gR on the rhs of 
the flow equation ([6]) generated by the Hessian of the 
Polyakov term in ([7]) . We report the details of the com- 
putation of the Hessian in the appendix. Following the 
notation of @, we find: 



+ C/y + 0{K^) , 
where V^ and Uij are as in ^ and 






V„V/3 

(9) 



W,f = Ka 



\5,,R+]-K'K^ +KX,K^^^ 



■yS 



^(^0 



(10) 



is the contribution to the Hessian coming form the 
Polyakov action. In ^ and PH)) we kept terms u" 
to order K^ and R. Following the procedure of 
to perform the functional trace and using the contrac- 
tion gapW"/ =Ka[2{D-3)R + 3K^] , finally gives the 
beta functions: 



dtfJ'k 



D 



Svr 



-Qi [GkdtRk 



dtKk = ^Q2[GldtRk]nk-'^Q2[GldtRk]KA 

D — 2 1 

dtfik = — — Qo [GkdtRk] - T-Qi [GldtRk] Mfc 



D-3 

327r 



Q2 [GldtRk] Ka . 



In pT|) the Q-functionals are defined as 

^h{z) n > 



^n[n\-<^ (-l)"/i(")(0) 



71 < 



(11) 



(12) 



and we introduced the fourth-order regularized propaga- 
tor: 



Gkiz) 



1 



KkZ^ + likZ + Rk{z) 



(13) 



The beta functions ([TT|) are valid for general embedding 
dimension D and for any admissible cutoff shape function 
Rk{z). They are the main result of this section. 

The effect of the inclusion of the hexatic stiffness 
through the Polyakov term in ([7]) is a modified renor- 
malization of the bending and Gaussian rigidities. For 
Ka = we reproduce the system of beta functions found 
in [a,13- Note that the additional renormalization of the 
Gaussian rigidity vanishes in the physically relevant case 
£) = 3. 



D. Crinkled phase 

We can start to consider the beta function for the bend- 
ing and Gaussian rigidities when /i^ — 0. Employing the 
cutoff shape function Rk{z) = (fc"* — z^)9{k'^ — z) in pT|) 
we find: 



dtHk 



in 



D 



3Ka\ 
4. Kk J 



OtKk = —Z 1- —7 KAttk ■ 



(14) 



6tt 167r 

The beta function of Kk agrees with [31j while the beta 
function of Rk is new. In terms of the couplings ak = I/k^ 
and ak = l/Sfe we find: 

dtak = -^(o-^KAak 

4:TT \ 4 

-2 

dtak = -^{D-8 + iD- 3)KAak) . (15) 

DTT 



We observe that the system ([T5|) has a non-trivial fixed- 
point at: 



AD 1 
~3~K^ 



a* = 0, 



(16) 



together with the Gaussian one a* = a* = 0. Note 
that a* depends inversely on Ka, thus the non-Gaussian 
fixed-point can be controlled in a perturbative expansion 
for large hexatic stiffness. 



J 



When we re-introduce the surface tension and use the 
dimensionless variable /i^ = fe^/zfe the flow in the (afc, jlk) 
plane is described by the following beta functions: 



dtjlk = 



dtak 





, ^ 2 + ,1,-^4 + A^,^^ 

Ing i 






-^ ^2 + A, + ^4 + A^,-^ 




ak {D - iKAak) 


2{1 - ak) + ak~^k 1 Afe ^^_^2 + ^fc-^4 + A^- 


4 


(17) 


2-(4 + A^.-^) 


[ '^'''~^' \/4 + M?.-4 ^2 + A, + ^4 + M^ 


4 



We can recover the beta function of the dimensionless 
surface tension given in [2j by keeping the terms of order 
a jlk in (fT7|) . The flow is depicted in Figure 2 in the case 
Ka = 1 and D = i where the non-Gaussian fixed-point 
has values a* = 4 and /x* — 0.107, and is characterized by 
an attractive and a repulsive direction. Thus the inclu- 
sion of the Polyakov term in the effective average action 
([7]) had the effect of creating a non-trivial fixed-point, 
giving an explicit example of how a non-local term can 
alter the phase portrait of a model. 

The fact that non-local terms in the action modify the 
RG flow of the local coupling, without inducing a self- 
renormalization, is reminiscent of the effect WZWN term 
has on the flow of the NLSM coupling fl6|. 

Hexatic membranes have a continuous phase transition 
in two dimensions and are characterized by a continu- 
ous symmetry (re-parametrization invariance), but since 
the phase transition is driven by long-range Coulomb in- 
teractions between Gaussian curvatures, there is no con- 
tradiction with the Mermin-Wagner-Hohenberg theorem 
|17| . The fixed-point value a, depends continuously on 
the hexatic stiffness Ka and one actually have a line of 
fixed-points as in the Kosterlitz-Thouless phase transi- 
tion Q. As we will see in a moment, also the critical 
exponents depend continuously on Ka- 

At the non-Gaussian fixed point the membrane has 
a non-trivial fractal dimension. This is related to the 
mass critical exponent by the general relation d_F = - [l| , 
where v is minus the inverse of the negative eigenvalue 
of the stability matrix. Linearizing the flow around the 
non-Gaussian fixed-point gives the estimate: 



dp = 2 



D~2 D 



TT 3Ka 
D~2 f D \i 



2{D_i^_^( D Y 



D 

3i^ 



(18) 



where /i* is the fixed point value of the dimensionless 
surface tension. The first correction term in (jlSp agrees 
with the one found in [2]. The fractal dimension is shown 




Figure 3: The fractal dimension dp as a function of 1/ Ka for 
D = ?,. The upper curve is the leading order term in 1/Ka 
while the lower curve is our estimate (I18II. 



in Figure 3 as a function of the hexatic stiffness for _D = 3; 
since its value is bigger than the classical dimension d = 2, 
the membrane is crinkled, i.e. it is infinitely rugged but 
still spatially extended. Our estimate (fT5|) is strictly lower 
than the leading 1/Ka result for all values of the hexatic 
stiffness and it tends to dp — 2.709 as Ka — > 0. 



E. Conclusions 

In this paper we studied hexatic membranes using a 
geometrical approach based on the effective average ac- 
tion. The 0(2)-model coupled to the fluid membrane 
has a non-trivial RG fixed-point where the membrane 
is crinkled and has fractal dimension depending on the 
hexatic stiffness and on the embedding dimension with 
values between two and three. 

This example shows how matter coupled to fluctuat- 
ing geometries can induce physically interesting phases 
through the generation of non-local invariants. The cal- 
culation of the present paper gives an example of the non- 
trivial role played by non-local terms in the effective ac- 
tion, showing that their influence can drastically change 
the phase structure of a given model. In this light, it is 
important to consider the influence of non-local terms in 



other models of fluctuating geometry, such as membranes 
in higher dimensions J18i] or quantum gravity |19| . 

This work also opens the road to the study of matter- 
membranes systems by the methods of the effective av- 
erage action and the relative exact flow equation. The 
question that naturally arises is how flat space universal- 
ity classes of 0(A'')-models [20], other than N — 2, are 
dressed by the geometrical fluctuations of the membrane. 
Are there infinitely many fixed-points in the A^ = 0, 1 
cases? Does the Mermin-Wagner-Hohenberg theorem ap- 
plies to the N > 2 cases? This, and the inclusion of topo- 
logical excitations in the N = 2 case, will be subject of 
further studies. 
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Appendix A 

We report here some details about the calculation of 
the contributions to the beta functions of k^ and Rk in- 



duced by the Polyakov action: 



lp[g]= J VaR^R- 



(Ai) 



Since we are expanding the rhs of the flow equation (|6]) 
to order K^ or R, it is enough to keep only contributions 
of order R in the Hessian of the Polyakov action (jAlf) . 
The expansion is conveniently written as: 



S^Ip[g]^ 2 I cfx^SR^SR 
+ 2 I d^Xy/gR^SR 



+ 0{R^,RK^), (A2) 



Sgfi„=H^ 



with h,,i, = —2v^Kl,„ the first variation of the metric 



ly'iyiK'^Ki!^^) the second. After 



and H^^ — 2{day^dpu'^ 

the expansion is performed, we can set the background 
metric, and implicitly the embedding field, equal to the 
metric of the two dimensional sphere; in this way we make 
all derivatives of the curvatures vanish, but we will still 
be able to disentangle the operators R and K^ . Using 
the variation 5R = g'^-'ASg^,^ + V^V^^^ - y'^'Sgf^^R 
we find: 



J 



Jd^Xy^R^SR 



d^Xy^SR—SR 



d^x^Hf,^iRg^"' + V>'V^R] +0{Ri\RK'^'' 
d^x^^ i^'Aiy' + 0{R'^,RK^) 



d^x^ f hAh + 2/iV^V'' V + /ia/S^V^V^V^V V 
d^^Vg (l^'^^ + KpK^"^) ^'At.^ + 0(i?2) . 

I 



0{R') 



(A3) 



In (IA3I) we made the substitutions: 



1. 



2 ^ 

8 



that are allowed, to this order, under the trace in the flow 
equation and we used the relation KapK°'^ = K^ — R 
to simplify the flrst term. Inserting (jASp in (jA2[) gives 
the contribution (|T0)) of the Polyakov action (jAip to the 
Hessian © of the effective average action. 
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